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ABSTRACT 

We investigate the spatially-resolved star formation relation using a galactic disk formed in a comprehensive 
high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation 
as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation 
rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, 
reducing the star formation rate. Because we have self-consistently calculated the location of ionized gas, 
we for the first time are able to make spatially-resolved mock observations of star formation tracers, such as 
Ha emission. We can also observe how stellar feedback manifests itself in the correlation between ionized 
and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3 x 10 M©, we find that the 
correlation between star formation rate density (estimated from mock Ha emission) and H2 density shows 
large scatter, especially at high resolutions of < 75 pc that are comparable to the size of giant molecular clouds 
(GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution. By 
examining the evolving environment around star clusters, we demonstrate that the breakdown of the traditional 
star formation laws of the Kennicutt-Schmidt type at small scales results from a combination of stars drifting 
from their birthplaces, and molecular clouds being dispersed via ionizing radiation and supernova feedback. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: star clusters — galaxies: structure 
— stars: formation — stars: evolution — ISM: structure — ISM: evolution 



1. INTRODUCTION 

Understanding how galaxies have evolved into their current 
shapes is one of the most intensely pursued, yet one of the 
most poorly understood topics of contemporary astrophysics. 
By definition the process of galaxy formation cannot be sepa- 
rated from the problem of star formation. Any analytic, semi- 
analytic, or numerical framework to build galaxies requires a 
model to form stars and describe their feedback. But, how 
good is our understanding of star cluster physics, especially 
when the star formation and feedback themselves cannot be 
understood outside the galactic context? 

In numerical studies, one may verify the star formation 
models by employing them in galactic simulations, and then 
comparing the resulting galaxies with the gold standard we 
already possess: observed galaxies. As a matter of fact, we al- 
ready have an abundant supply of galactic observations accu- 
mulated over decades to which we can compare our simulated 
galaxies, including photometric images of galaxies at various 
wavelengths, spectral energy distribu tions, star forma tion his- 
tories, star formation relation s (e.g. lKennicuttlll998l: Kenni- 
cutt et al. 2007), and so on. Yet, in order to make a reliable 
and meaningful comparison between them for the purpose of 
understanding galactic star formation, one must ensure that 
the galaxy simulation be properly coupled with star cluster 
physics. To that end, the following conditions should be met: 
(a) Sufficient resolution in time and space to describe the rel- 
evant star cluster physics. Resolution of < 10 pc is needed to 
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resolve the scale of star-forming molecular clouds and to pro- 
duce a multiphase interstellar medium (ISM). A correspond- 
ing timestep would res olve the lifetime of massive stars and 
molec ular clouds (e.g. iTevssier et all 120 1 Ol : fBournaud et alj 
2010). (b) Physically-motivated models for star cluster for- 
mation and feedback with as little subgrid physics as possible. 
Models should correspond with the spatial and temporal reso- 
lution of the simulation. For example, a star formation model 
that solely depends upon the star f ormation relatio n observed 
at ~ kpc scales (e.g. Esfr ~ £g a 4 ; lKennicutlll998h may not 
be desired in high-resolution simulations. Similarly, a phe- 
nomenological feedback model such as stopping the cooling 
near a star cluster particle right after its birth might not be a 
suitable choice given fine temporal resol ution (see, however, 
Hopkins e t alJl20lH : IStinson et alJl2012h . Not surprisingly, 
simulations at this level of sophistication require enormous 
computational costs. For this reason one may arguably state 
that no simulation has reliably checked the physics of galactic 
star formation against observed galaxies. Indeed, galactic star 
formation still remains to be fully understood. 

However, thanks to the expeditious developments of both 
numerical techniques and computer hardware, simulations 
with significantly improved resolution are now feasible at in- 
creasingly cheaper costs. Utilizing adaptive mesh refinement 
techniques, in the first paper of this series dKim et alj|2012l 
hereafter Paper I) we simulated the behavior of star clusters 
and their embedding galaxy in a single self-consistent nu- 
merical framework. The implementation includes a new way 
of describing feedback from star-forming molecular cloud 
(SFMC) particles by combining ionizing stellar radiation and 
supernova explosions. The radiation feedback from each of 
< 10000 SFMC particles is rendered by tracing the ultraviolet 
photon rays on the fly. This means that our approach does not 
rely on a subgrid model to handle the stellar radiation feed- 
back, but actually solves the transfer equation. Joined with 3.8 
pc resolution, this new scheme allowed the authors to study 
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TABLE 1 
Simulation Suite Description 



Physics 2 


MC-TF 


MC-RTF 


Star-forming molecular could (SFMC) particle formation 


o 


o 


SFMC feedback: supernova explosion 


o 


o 


SFMC feedback: ionizing radiation 


X 


o 



"For detailed explanation, see S|2]or Paper I, o = included, x = not included. 

the escape of ionizing photons from star-forming clumps and 
from a galaxy. By simulating a galactic halo of 2.3 x lO n M , 
we found that the galactic escape fraction is dominated by 
a small number of molecular clouds with exceptionally high 
escape fraction. We also discovered that the escape fraction 
from a SFMC particle rises on average from 0.27% at its birth 
to 2.1% at the end of a cloud lifetime, 6 Myrs. 

The machinery established in Paper I is very promising and 
provides us with an one-of-a-kind look into the physics of star 
formation and its intertwined nature with the galactic context. 
Moreover, our realistic stellar feedback model enables us to 
make unprecedented simulated observations of star formation 
tracers, such as Ha emission. In the second paper of the se- 
ries, we use the numerical framework developed in Paper I 
to carry out an analysis of self-regulated star formation in a 
dwarf-sized galaxy. We make a direct comparison of star for- 
mation tracers between simulated galaxies and observed ones. 
Special attention is paid to the spatially-resolved star forma- 
tion relation observed in local galaxies, and why it breaks 
down at small scales. We demonstrate how galactic star for- 
mation can be better understood by significantly expanding 
the usage of cutting-edge numerical simulations. 

This article is organized as follows. The physics in the sim- 
ulation code and the initial condition are briefly explained in 
Sj2]and $3] respectively. ^overviews me global star formation 
of the performed runs. fJ3]is devoted to the spatially-resolved 
star formation relation via simulated maps of physical observ- 
ables such as ionized gas tracers, Ha. We also inspect the 
evolving environment of star-forming gas clumps. Assembled 
in fj6]are the summary and conclusions. 

2. PHYSICS IN THE CODE 

Our enhanced version of the Eulerian adaptive mesh re- 
finement Enzo-2. 1 5 dBrvan & Normanlll997l iNorman et al.l 
120071: 1 Wise & Abelll2011b iKim et al.1 120091120111) contains 
all relevant physics previously considered in galactic simula- 
tions as well as the new treatments of SFMC particle feed- 
back discussed in detail in Paper I. The framework developed 
in Paper I is designed to investigate galactic star formation of 
a dwarf-sized galaxy in a self-consistent simulation, without 
omitting the detailed interaction between the SFMC particles 
and their surrounding galactic context. While the key ingre- 
dients of the framework are carefully described in Paper I, we 
briefly review its major features focusing on the SFMC feed- 
back model. Table [TJ also summarizes the feedback models 
we have employed in the reported study. 

The ZEUS hydrodynamics module in Enzo is employed to 
solve the Euler equations for collisional fluid (Stone & Nor- 
man Il992allbh . We use Enzo's multi-species non-equilibrium 
chemistry and cooling modules to track reactions among six 
species (H, H + , He, He + , He ++ , e~) and radiative losses. 
Cooling by metals is included in gas above 10 4 K using the 

5 http://enzo-project.org/ 



cooling rates tabulated by Sutherla nd & Dopltal (119931) . and 
in gas below 10 4 K us ing th e cooling rate approximation 
of iKovama & Inutsukal (12002b . The mass thresholds above 
which a cell is refined by factors of 2 in each axis are functions 
of a refinement level / as M' mf gas = 2-°- sm ■ 4.87 x 1O 6 M 

for gas, and M' ef part = 2-°' 533/ • 4.87 x 1O 6 M for particles, 
respectively. At the finest level / = 13 of Ax = 3.8 pc reso- 
lution, a cell would have been refined if more than 3000 M & 
in gas or 40000 M ~ 2.5 Mdm in particles. This resolution 
is in accord with the Jeans length for a dense gas clump of 
n = 1000 cm~ 3 at ~ 100 K, at which point a corresponding 
Jeans mass of 2000 M collapses to spawn a SFMC particle. 
Indeed a finest cell produces a SFMC particle of initial mass 
Mjjjc = 0.5p gas Ax 3 when (a) the proton number density ex- 
ceeds the threshold 1000 cm -3 , (b) the velocity flow is con- 
verging, (c) the cooling time f coo i is shorter than the gas dy- 
namical time fdyn of the cell, and (d) the particle produced has 
at least M t h res = 1OOOM . This formulation guarantees that 
a SFMC particle is created before an unresolved gas clump 
interferes with the consistency of hydrodynamics. 

In our simulation, the later evolution of a SFMC particle 
describes three key aspects of star formation in Giant Molec- 
ular Clouds (GMCs): (a) 0.24 is instantly turned into 
stars, but the rest is considered as atomic or molecular gas 
that does not participate in star formation, (b) Among stars, 
0.04 MfJJc lS massive enough to commence Type II supernova 
pfeaking at the particle age of 5 Myr. The rest of the stellar 
mass, 0.2 M™^, is considered as long-lived stellar mass, M*. 
(c) 0.76 + 0.04 = 0.8 Mjjjj-. is eventually recycled back into the 
ISM with the energy of exploding supernovae (7.5 x 10~ 7 of 
the rest mass energy of 0.2MJJ}c). 2% of the ejected mass is 
considered as metals enriching the ISM. Overall, the mass of 
a SFMC particle evolves as 



M M c(t) 



M mit 
M MC 



1-0.8 [o.5-0.5erf{V8(5-r)} 



where erf () is the Gauss error function, and T = t — f cr is the 
particle age in Myr with particle creation time f cr . 

In order to describe the ionizing ultraviolet photons from 
young star clusters, we perform a self-consistent transport cal- 
culati on of monochromatic stellar radiation from SFMC par- 
ticles (lAbel & WandeltlHOQllWise & AbellloTl IKim et alJ 
1201 ll) . From T = to 6 Myrs after the birth, ionizing radi- 
ation luminosities Lmc(M) = ?Mc£phA^Mc(', t) are assigned 
to the SFMC particles that are more than 2.5 kpc away from 
the galactic center. Here ^mc = 6.3 x 10 46 photons s^'Mq 1 
is the lifeti me-averaged ionizing lu minosity per solar mass 
in clusters (Murray & Rahman! 1201 Oh . a nd £ pn = 16.0eV is 
the mean energy per deposited photon (Wh alen & Norman] 
I2006D . 6 The SFMC particles in the inner 2.5 kpc disk do 

6 By using the particle mass, Myic(Lt), rather than the stellar mass inside 
we aim to approximately encapsulate various other types of feedback beyond 
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FIG. 1 . — Left: face-on density-weighted projection of density in the MC-RTF run in a 20 kpc box at 13.3 Myr into the high-resolution evolution. High-density 
clumps which bear a number of SFMC particles are noticeable. Right: density-weighted projection of photoheating rates. Note that the photoheating rates are 
low in the galactic center (< 2.5 kpc in radius, dotted circles) because the star-forming molecular cloud (SFMC) particles there are not ray-emitting. 



not emit radiation so that we concentrate on the evolution of 
clouds in the galactic spiral arms and outer disk, and expe- 
dite the radiative transfer calculation. Each of the 12 x 4 3 
rays isotropically cast from the particle is traced until most 
of its photons are absorbed or until it reaches the edge of the 
computational domain, while being adaptively split into child 
rays whenever the angular resolution associated with it grows 
larger than 0.2 (Ax) 2 of a local cell. 7 Photons in the emitted 
ray interact with hydrogen in the surrounding gas in three dif- 
ferent ways: they (1) photoionize hydrogen (2) photoheat the 
gas with the excess energy above the ionization threshold, and 
(3) exert outward momentum when they are absorbed by the 
gas cell. Dust absorption, or the infrared photons re-emitted 
by UV-irradiated dust grains and the momentum imparted by 
them are not explicitly included in our calculation. 

3. INITIAL CONDITIONS 

We construct a dark matter halo of 2.3 x 10 ll M ( r ) that fol- 
lows the Navarro-Frenk-White profile with a concentration 
c = 10 and a spin parameter X = 0.05. Gas grids are generated 
by splitting particles with an initial baryonic mass fraction of 
10% and metallicity of Z; n i t = 0.003 Z . We then produce a 
galaxy embedded in a dark matter halo by performing a coarse 
resolution (15.2 pc) simulation for ~ 1 Gyr in a 32 3 root grid 
box of 1 Mpc 3 . When a relaxed, well-defined galactic disk has 
emerged, ~ 1.3 x 10 10 M Q is in 3.4 x 10 6 star cluster particles 
formed at low resolution, and ~ 6.5 x 10 9 M e is in gaseous 
form, all embedded in a halo of 1 .3 x 10 7 dark matter particles 
(16000 Mq each). We now employ high-resolution refinement 
criteria (3.8 pc, 13 additional levels of refinement) to resolve 
the galaxy down to the size of a molecular cloud, and make 
the entire galaxy evolve for ~ 30 Myrs. As a consequence, the 
computational expense of performing a high resolution calcu- 
lation for a galactic dynamical time is saved, but one may still 

photoionization, such as protostellar outflows, stellar winds, and radiation 
pressure. More discussion on this in Paper I. 

7 Other techniques to speed up the radiation calculation include: (a) A ray 
is no longer split when it is more than 5 kpc away from the source, a distance 
far enough for a dwarf-sized galaxy, (b) Two rays merge to form a single 
ray when the distance from the source is more than 10 times the separation 
between the two sources, (c) Two radiating SFMC particles merge if they are 
separated by less than 5 pc in space, and 3 x 10 4 yr in age. 



observe the star-forming gas clumps for their typical lifetime. 
To exclude the period in which the galaxy is trying to reach 
new equilibrium with freshly imposed high resolution, we in- 
vestigate the galaxy from 10 to 30 Myr into the evolution in 
the remainder of this article. Figure [T] shows the snapshot of 
our simulated galaxy in the MC-RTF run at 13.3 Myr into the 
high-resolution evolution. 

4. RESULTS I: GLOBAL STAR FORMATION 

We analyze two simulations listed in Table [T] to examine 
the global star formation, and the spatially-resolved star for- 
mation relation in a simulated dwarf-sized galaxy. First we 
briefly overview the global star formation rate on the galactic 
disk and the structure of the galactic ISM. 

4. 1 . Self-regulated Galactic Star Formation 

Shown in Figure [2] is the instantaneous galactic star for- 
mation rate, dM^/dt (with dt = 0.37 Myr), of the run with 
only supernova feedback, MC-TF, and the run with both ra- 
diation and supernova feedback, MC-RTF. During the ~ 20 
Myrs evolution with 3.8 pc resolution that is analyzed, the 
star formation rate (SFR) on average is 4.6OM yr _1 and 
3.57 M Q yr , respectively. In other words, the MC-RTF run 
forms 22.4% less stellar mass than the MC-TF run does. Our 
realistic description of stellar feedback helps to self-regulate 
galactic star formation, and achieve reasonable star formation 
rates during the analyzed period. 

Figure [3] illustrates the relationship between the time- 
averaged SFR surface density, Esfr, an d 8 as surface density, 
£ gas , for the two runs between 13.3 Myr and 30.7 Myr of the 
high-resolution evolution. Each of the densities is measured 
on a galactic disk of 9 kpc in radius excluding the central 2.5 
kpc in which SFMC particles do not radiate ionizing photons 
(see $2). We estimate the SFR surface density by using the 
stars of age less than T = 0.37 Myr. The movement of the 
datapoints in time on this plane is relatively small during the 
analyzed period which covers only 17.4 Myrs. The error bars 
here show the extent of such movements; and the movement 
in x-axis (variation in E gas ) is so small that its error bar is not 
visible in the plot. We acknowledge that the overall suppres- 
sion by stellar radiation is not large enough to move the dat- 
apoints down to the observed global Kennicutt-Schmidt rela- 
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FIG. 2. — Time evolution of the instantaneous galactic star formation rate in 
the unit of M®yr from 13.3 Myr to 30.7 Myr of the high-resolution evolu- 
tion. For the MC-TF (solid line, only with supernova feedback) and MC-RTF 
runs (dashed line, with both supernova and stellar radiation feedback). 

tion dKennicuttl lT99 8 ) . The suppression ratio could have been 
larger if the SFMC particles in the inner disk had also radi- 
ated ionizing photons (< 2.5 kpc from the galactic center) in 
which the largest fraction of stars tend to spawn in general. 
We will come back this issue in f|6] 

4.2. Interstellar Medium (ISM) Structure 

We now focus on how ionizing radiation feedback regulates 
star formation. Photoionizing radiation from massive stars is 
one of the important factors that drives changes in the environ- 
ment of star-form i ng clumps, and self-regulates star formation 
dWhitworthl 1 9791 : lMatznerl2002l) . In order to characterize the 
nature of the impact of ionizing stellar radiation feedback, we 
compare the galactic ISM structures of the two runs, MC-TF 
and MC-RTF. The bottom rows of Figure [4] displays the two 
dimensional probability distribution functions (PDFs) of pro- 
ton number density and gas temperature, colored by gas mass 
in each bin. The measurement is made in the sphere of 10 kpc 
radius centered on a galactic center (but excluding the inner 
2.5 kpc sphere) at 13.3 Myrs into the high-resolution evolu- 
tion for each of the performed simulations. 

The ionizing radiation from young SFMC particles heats 
the surrounding dense gas up to ~ 10 5 K and pushes the PDF 
in an upward direction. In the MC-TF run dense gas cells 
with number density above 10 3 cm~ 3 , the threshold density 
for the SFMC particle formation (see Sj2), are almost com- 
pletely turned into SFMC particles. However, in the MC-RTF 
run, zone "A" still hosts very dense gas cells above 10 3 cm~ 3 , 
which are beyond the threshold density. Photoheating by stel- 
lar radiation retains such hot dense gas which otherwise would 
have been unstable against Jeans fragmentation and deposited 
into SFMC particles. At the same time, zone "B" depicts rel- 
atively less dense gas cells heated up to ~ 10 5 K, but mostly 
clustered around ~2x 10 4 K just above the bump in the cool- 
ing curve by the [OIII] forbidden line. The top rows of the 
same figure show one dimensional PDFs of gas number den- 
sity at the same epoch, binned by cold (< 2 x 10 3 K), warm 
(2 x 10 3 - 2 x 10 s K), and hot gas (> 2 x 10 5 K). In the plot 
for the MC-RTF run, the increased masses of the warm gas 
above 10 3 cm~ 3 (corresponding to zone "A") and the hot gas 
at around ~ 10~ 2 cm~ 3 are prominent. 



FIG. 3. — Time-averaged star formation rate (SFR) surface density, Zsfr. 
and gas surface density, £ ea s, between 13.3 Myr and 30.7 Myr of the high- 
resolution evolution, for the MC-TF (red circle) and MC-RTF runs (blue 
square). Also plotted are the observed global Kennicutt- Schmidt relations: 
the best fits to the disk-averaged galactic observations by Kennicutt 1 1998, 
long dashed line w ith a shaded region representing their errors) and by 
IBouche et al. (2007, dot dashed line that includes high-z galaxies). 

5. RESULTS II: SPATIALLY-RESOLVED STAR 
FORMATION RELATION 

In this section we move to the spatially-resolved star for- 
mation relation on a simulated galactic disk. We start by de- 
scribing our methodology which for the first time allows us 
to make spatially-resolved mock observations of galactic star 
formation tracers, such as Ha emission. 

5. 1 . Simulated Observation of Star Formation Relation 

The traditional star formation tracer of Ha emission 
(B aimer line of n = 3 2, 1.88 eV) maps out the ionized 
gas around young massive stars. The site of highly ionized 
gas, however, may not overlap with the site of dense collaps- 
ing gas in which on-going star formation takes place, usually 
identified by peaks of molecular hydrogen, H2. This discrep- 
ancy may become even more problematic when the observed 
galaxy is resolved with higher spatial accuracy which distin- 
guishes such regions from one another. Therefore in an ob- 
servation with high spatial resolution, one may question the 
validity of the traditiona l star formation tracers, in particu- 
lar, Ha emission (e.g. lOnodera et "all 1201 Ol : ISchruba et alj 
l2010t iKomugi et al.ll2012l) . To understand and characterize 
this phenomenon, we now produce mock observations of var- 
ious physical properties in our simulated galaxy. Our galaxy 
is an ideal laboratory for such an exercise not only because 
we have self-consistently tracked radiation fields like pho- 
toionization rates prompted by < 10000 SFMC particles with 
non-equilibrium primordial chemistry, but because the galaxy 
is simulated with high resolution, 3.8 pc, better than any re- 
ported observation of local galaxies. A high dynamic range in 
our simulation also facilitates improved estimation of molec- 
ular hydrogen density. In the remainder of f|5] we mainly fo- 
cus on the MC-RTF run in which radiation fields are tracked 
through transport calculation, unless specified otherwise. 

To begin with, we describe how the SFR surface density 
is mapped from a mock observation of Ha emission. In the 
hot HII regions around young massive stars, the main process 
contributing to the Balmer lines is the recombination of ion- 
ized hydrogen, H + . Therefore, first and foremost, we estimate 
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FIG. 4. — Bottom rows: joint probability distribution functions (PDFs) of gas number density and temperature colored by mass in each bin, for a 10 kpc sphere 
from a galactic center (but excluding the inner 2.5 kpc sphere) at 13.3 Myrs into the high-resolution evolution (left: MC-TF run, only with supernova feedback, 
right: MC-RTF run, with supernova and radiation feedback). The vertical dashed line in each plot denotes the density threshold for SFMC particle formation. 
The ionizing photons from SFMC particles heat the surrounding dense gas up to ~ 10 5 K and prevent star formation. It increases the amount of very dense gas 
which otherwise could have been gravitationally unstable (> 10 3 cm -3 ; zone "A"). Top rows: one dimensional PDFs of gas number density that share the x-axes 
with the bottom rows. In the MC-RTF run, a significant amount of warm gas (2 x 10 3 — 2 x 10 5 K) is existent above 10 3 cm -3 . 



the Ha emission by recombination per volume element dV as 

dL Ho , tR =47tj lia ,Rd v ( 2 ) 

/i^-wo el „ | e\ — 26 rri — 0.942 — 0.031 InT/t 
= 47TX 2.82x10 T 4 *n e n H +dV{3) 

in erg s , where jua,R refers to the Hot emis sion rate by re- 
combination (Eq.(9) of Don g & Draindl201 ll) . n e and « H + 
to the number density of e and H + , respectively, and T4 to 
the temperature T e in the unit of 10 4 K. Second, we add the 
enhancement of Ha emission by collisional excitation. This 
boost is usually minor because, in typical HII regions, the 
neutral hydrogen fraction is very small, and few free elec- 
trons have sufficient ener gy to excite bound el ectrons of neu- 
tral hydrogen to n > 3 (Lurid iana et al.ll2003h . We estimate 
the contribution of collision excitation n = 1 — > 3 followed by 
de-excitation n = 3 — > 2 as 

dLn at c=E 32 qi3n e njidV (4) 



^0l e - 12 - uv ^n e n H dV (6) 



= 1.30 x 10~ 17 



where Eij is the energy associated with the tran- 

h <7b 



in erg s 

sition ; — > j, qu is the rate coefficient of excitation n = 
1 — > 3, ri3(r e ) is the Maxwellian-averaged effective collision 
strength, 8 ©i = 2 is the statistical weight of the originating 

8 r 13 (r c ) =0.350- 2.62 x 10- 7 r c -8.15 x lO" 11 r c 2 + 6.19 x io- i5 t c 3 



level, Icq is the Boltzmann constant, m e is the mass of an elec- 
tron, and (Inti'/knmiyi 2 = 8.63 x 10~ 6 . Here we fix T e at 
8000 K, a characteristic temperature inside a typical HII re- 
gion that is poorly resolved in our simulation. 9 Third, we 
subtract the emission from the diffuse ionized gas (DIG) by 
removing the contribution by cells whose proton number den- 
sities are smaller than 10 cm -3 or whose temperatures higher 
than 10 5 5 K dSchruba et alj|2010h . Lastly, we convert the lu- 
minosities in Ha emission into the star formation rates with 



dM*/dt (M Q yr" 1 ) = 5.3 x 10~ 42 L Ha (erg s 



-In 



(7) 



in which a correction for dust extinction is not considered (e.g. 
ICalzetti et "ail 120071) . In order to distinguish the SFR den- 
sities measured in these mock observations from the actual 
values in simulations (e.g. as in ij4.ll ). we hereafter define 
PSFR,est = 5.3 x 10~ 42 dLua/dV and IsFR,est as the SFR (sur- 
face) density estimated from a simulated observation. 
We now discuss how the H2 surface density is mapped. 



for 4000 K < T e < 25000 K from Table 3 of Ag garwal fT983l) . 

9 Collisional excitation may falsely dominate Ha luminosities if the tem- 
perature is inadequately evaluated because Eq.{6) has a steep temperature 
dependence. In fact, individual HII regions in the reported simulation are not 
properly resolved, and monochromatic photons (16.0 eV) from SFMC par- 
ticles often keep overly-heated temperature (> 2 x 10 4 K) in the neighbor- 
ing under-resolved HII regions. When this temperature is naively adopted, 
Eq.{6} might overestimate Lhh.c- I n an attempt to compensate for our poor 
resolution, the temperature in Eq.j6} is thus fixed at 8000 K, a characteristic 
temperature inside a typical HII region. With this correction, the contribution 
by collisional excitation is indeed found to be only marginally important. 
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FIG. 5. — Left: SFR surface density, £sFR.est, in the unit of W. yr 1 kpc 2 evaluated from a mock observation of Ha emission using Eq.jjJ in the MC-RTF 
run. Shown at 13.3 Myr into the high-resolution evolution in a 3 kpc box centered on a clump that actively forms SFMC particles. These face-on projections are 
generated with uniform 15 pc resolution. Locations of 20 local maxima are marked with black circles. Right: H 2 surface density, £H 2 ,est, in the unit of M e pc~ 2 
evaluated with Eq.(8). Again, locations of 20 local maxima are marked with circles. Note that Ha and H2 peaks are closely related, but do not perfectly coincide. 



Molecular hydrogen is not explicitly traced in our chemistry 
solver which follows only six atomic species. We instead use 
an approximation to estimate the molecular hydrogen frac- 
tion, fn 2 , as a function of gas density and metallicity by 



fn 2 



I- 3 - 



4 l+0.25s 



with parameters defined as 



ln(l +0.6* + 0.01 J 2 ) 



/ 0.365 1 



0.6 T c 

j~0.77(l+3.1Z 



T c = ~ 321 Z - 

Mh Vp gas 



Pgas 



(8) 

(9) 
(10) 

(11) 



where approximations in Kuh len et alj d2012h have been ap- 
plied to the original formula of the scaled radiation fie ld % 



(IKrumholz et"dl200ai2009tlMcKee & Krumholzl20ioh . and 
an approximated column density of the cell, E = p 2 /| Vp |, has 
been employed for the cell' s dust optical depth T c (the S obolev 
approximation described in IKrumholz & Gnedinll201 lb . Here 
Z' is the metallicity normalized to the solar value, Od is the 
mean extinction cross section by dust per hydrogen nucleus, 
jIIh is the mean mass per hydrogen nucleus, and p 2as is in 
the unit of g cmT 3 . We hereafter define pH 2 ,est = /iuPgas an d 
^-H 2 ,est as the H2 (surface) density evaluated from a simulated 
galaxy. This is to emphasize that the our H2 density is an 
approximate estimate, not a real value which is explicitly fol- 
lowed in simulations. 

Because the location of ionized gas is self-consistently 
computed in our high-resolution radiation hydrodynamics 
simulation, we are able to investigate the correlation between 
ionized and molecular gas in exceptional detail. With psFR.est 
and Pn,.est estimated as above, projections of both densities 
can be readily generated. One such example is shown in Fig- 
ure[5]where a face-on projection of SFR density, EsFR,est7 an< ^ 
that of H2 density, £H 2 ,est> are produced and compared side 




0.5 1.0 1.5 2.0 2.5 

!og[ Approx. H 2 Surface Density (MsurVpc~2) ] 

FIG. 6. — The distribution of local maxima found in Figure|5]on the plane 
of SFR surface density, ZsFR.est (estimated from mock Ha emission), and 
approximate H 2 surface density, EH 2 ,est- Peaks of EsFR.est are in blue and 
peaks of Ej^.est in red. A blue star represents the average ZsFR.est Cv-axis) 
and the average EH 2 ,est (jc-axis) for the peaks in the ZsFR.cst map, while a 
red filled c ircle is for the ZHi.es t peaks. These markers are chosen to mimic 
Figure 4 of Schruba et al. 1 2010). Plotted as a dashed line is the best fit to the 
spatially-resolved observation of M51 byfKennicutt et al. 1 20071). 

by side. These maps are generated with uniform 15 pc reso- 
lution 10 at 13.3 Myr into the high-resolution evolution of the 
MC-RTF run, in a 3 kpc box centered on a gas clump that 
harbors multiple young SFMC particles. Also displayed are 
the locations of 20 local maxima of each map marked with 
black circles. 11 The resolution of these maps can be regarded 
as a mock observation equivalent of an aperture size (diame- 
ter) in photometric observations of local galaxies. For Figure 
[5] the resolution of 15 pc is chosen in order to demonstrate 
how the locations of peaks in one map differ from the ones in 
the other. In these images, it is unmistakable that Ha peaks 



10 Using FixedResolutionBuf f er in yt I Turk et al. 201 1 ) that deposits 
meshes of varying sizes into a pixelized, fixed-resolution array. 

1 1 Using maximum_f ilter in scipy package at http://www.scipy.org/. 
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FIG. 7. — Estimated SFR surface density (left, in Mq yr~' kpc -2 , from mock Ha emission) and H2 surface density (right, in M . pc~ 2 ). Shown at 13.3 Myr 
into the high-resolution evolution of the MC-RTF run in a 20 kpc box centered on the galactic center. The same mock observation te chni que as in Figure [5] is 
employed with uniform 50 pc resolution. The central region of radius 2.5 kpc is cut out in these maps (inner dashed circle; see [|2]and 35.2j . The locations of 50 
local maxima are marked with black circles. Note that star-forming clumps are highly clustered; compare with the projected density in the left panel of FigurefT] 



and H2 peaks are closely related, but they do not always co- 
incide perfectly. This experiment verifies the proposition we 
discussed in the beginning of this section. That is, the site of 
highly ionized gas (typically identified with Ha peaks) may 
not necessarily overlap with the site of on-going star forma- 
tion (typically identified with H2 peaks). 

To portray this discrepancy in a clearer way, in Figure|6]we 
place the local maxima of Figure [5] on the plane of SFR sur- 
face density, EsFR.esb an d H2 surface density, £h 2 est- Both 
axes represent estimated values from the simulated galaxy, 
MC-RTF, at 13.3 Myr into the high-resolution evolution. Each 
data point can be regarded as a 15 pc aperture average at the 
local maxima in the map of either EsFR.est orE^est- SFR sur- 
face density peaks are marked with blue circles and H2 peaks 
with red circles. It is obvious that the EsFR.est peaks are not 
necessarily the peaks in Eldest ma P> ar, d y i ce versa. Further, 
some data points suggest that EsFR.est P ea ks could have very 
low £h 2 est ( u PP er l e ft m the P a nel), and vice versa (lower 
right in the panel). In other words, the locations of highly ion- 
ized gas do not precisely correlate with those of dense molec- 
ular gas. The discrepancy can be encapsulated by a separation 
between the averages of the two sets of peaks. A blue star 
represents the average EsFR.est an d ^H^.est f° r the peaks of the 
Esfr. est map, whereas a red filled circle represents the average 
for the peaks of the Ejj 2 est map. The separation between the 
blue star and the red filled circle characterizes how well Ha 
performs as star formation tracers. The separation becomes 
particularly large at high resolution, i.e. small aperture size. 
We will further discuss this issue later in £15.21 

5.2. Spatially-resolved Star Formation Relation 

We now apply the mock observation techniques described 
in <2j]to the entire galactic disk. With these techniques, it 
is possible to directly compare star formation tracers between 
simulated galaxies and observed ones. We can also observe 
how stellar feedback manifests itself in the correlation be- 
tween ionized and molecular gas at various resolutions. 

Figure|7]shows the same estimation of SFR surface density 
and H2 surface density, as in Figure [5] Instead of focusing 



on a single star-forming clump, Figure [7] looks at the galactic 
disk of radius 9 kpc (outer dashed circle) that encompasses 
most of the active SFMC particles. The uniform resolution of 
these maps, 50 pc, can be viewed as a diameter of apertures in 
our simulated observation. As in Figure [T] the central region 
of radius 2.5 kpc from the galactic center is specified with a 
small dashed circle. We exclude the contribution by this re- 
gion because the SFMC particles there have not emitted ioniz- 
ing photons by design (see Sj2]i and may thus contaminate our 
analysis. The locations of 50 local maxima are marked with 
black circles, which tend to cluster in a few small areas. This 
mock observation carefully mimics the highest resolution ob- 
servati on of a local galaxy, M33, performed by Schru baet alJ 
(2010, compare with their Figure 1). 

From our high-resolution galaxy data, we can make galactic 
maps like Figure|7]with arbitrary resolution (aperture size) as 
long as it is larger than the simulation's spatial resolution, 3.8 
pc. In the subsequent analysis, we analyze such maps gener- 
ated at seven different resolutions: 1200, 600, 300, 150, 75, 
50, and 30 pc. We then distribute all the pixelized elements 
of each galaxy map on the plane of SFR surface density and 
H2 surface density. The plots resulting from this process are 
in Figure [H] where we show how the pixels on a galaxy map 
are spread with respect to their SFR surface densities and H2 
surface densities. The panels represent four different resolu- 
tions of the map, and each of the data points is colored by 
density-weighted projected temperature. Additionally, data 
points with highlighted edges stand for up to 50 local peaks 
in EsFR.est (points with blue edges) and EH 2 , es t (points with 
red edges). 12 A blue star and a red filled circle each repre- 
sents the average location of each set of peaks. 

12 ISchruba et S\ jlOlOh first identified CO regions and Ha peaks using 
their smallest aperture, d = 75 pc. Then an aperture of diameter [1200, 600, 
300, 150, 75 pc] is centered on each CO and Ha peak to measure the fluxes 
within that aperture. In order to reproduce such techniques in a simulated 
post-process, we first identify 50 peaks at 50 pc resolution. Then for a series 
of resolutions [1200, 600, 300, 150, 75, 30 pc], we find up to 50 local maxima 
that are nearest to the peak locations at 50 pc resolution. We use these peaks 
in the remainder of the analysis (e.g. in the making of Figure[8]and|9j. 
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FIG. 8. — The distribution of pixels of galaxy maps (e.g. Figure|7J on the plane of SFR surface density (from mock Ha emission) and H9 surface density. Each 
panel renders 1200, 300, 75, and 30 pc aperture averages at 13.3 Myr into the high-resolution evolution. Points are colored by projected temperature weighted 
by density. Also displayed are local peaks in ZsFR.est (with blue edge) and Eh-.. est (with red edge). A blue star and a red filled circ le each represen t s the a verage 
location of each set of peaks. A red dashed line shows a best fit in each panel, while a black line denotes a best observational fit by Kennicutt et al. 1 20Q3). 



When compared with observations, Figure [8] provides us 
with an illuminating view of the galactic star formation. First 
and foremost, one of the most unambiguous features in these 
panels is the large scatter in the data, up to ±3 orders of 
magnitude away from the best fit found in M5 1 using a 520 
pc aperture (black dashed line; iKennicutt et al.l l2007). espe- 
cially in the mock observations with high resolutions, < 75 
pc. The star formation relation observed on ~ kpc scales seem 
to rapidly break down once aperture sizes reach < 75 pc that 
is comparable to the scales of GMCs. This result is consistent 
with the recent observations of local galaxies b y many authors 
using smaller aperture sizes down to 75 pc (e.g. lOnoderaet al.1 
201Q: ISchruba"etini2l)T0l:lLiu et al.ll201lh . 

Each data point in these plots averages over star-forming 
regions with a r ange of evolutionar y stages and diverse Ha- 
to-H2 ratios (e.g. lSchruba et al.ll2.Q10h . Therefore, a data poin t 
with a large aperture (> 300 pc; e.g. Ke nnicutt et alj|2007h 
or a disk-averaged observation (e.g. lKennicutl fl998) yields 
not only a spatially-averaged estimate, but also a temporally- 
averaged one. In contrast, a data point with a small aperture 
(< 75 pc) captures only particular stages of GMC evolution, 
because the aperture is too small to average over a wide range 
of the evolutionary stages of GMCs. Hence, the correlation 
between ionized and molecular gas becomes loose at small 
scales. To further illustrate this effect, we paint various evo- 
lutionary stages of each data point (star-forming region) in 
Figure|8]by its projected temperature. In particular, in the 30 



pc resolution panel, a mild but clear transition is pronounced 
from lower right to upper left. It indicates a wide spectrum 
of star-forming regions on a single galaxy, from a cold dense 
region that actively forms stars (lower right) to a hot diffuse 
region that is likely dominated by the radiation from neigh- 
boring star clusters (upper left; hence exhibiting less active 
star formation). This supports the idea that at a size compa- 
rable to the scales of GMCs, each aperture re presents only a 
small fraction of evolutionary stages of GMCs ( Onodera et al.1 
l20ToLlSchrubaetal.l 2010; Ko mugi et al.ll2012h . We will re- 
investigate this issue in ^5.31 

Finally, we focus on the locations and averages of EsFR.est 
and EH 2 .est peaks, using the same technique and notation de- 
scribed in Figure [6] The bottom panel of Figure [9] shows 
the average locations of up to 50 local maxima identified at 
seven different resolutions. It is th e simulation equivalen t of 
the similar analysis performed by ISchruba et al.l d2010L top 
panel). As was examined in £ 15. H and Figure|6] the separation 
between the blue star and red circle encapsulates the discrep- 
ancy between the highly ionized gas and dense molecular gas 
at a given resolution. Although the locations of data points 
are not identical across the panels, 13 the separation between 
the two averages grows in both panels with smaller aperture 

13 This could be attributed to differences between M33 and our simulated 
galaxy, including the difference in evolutionary stage and mass. It could also 
be due to the incomplete ways to numerically estimate Ha and /h 2 . 
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FIG. 9. — Top: blue stars show the median Esfr and £h 2 f° r apertures centered on Ha peaks of M33. Red fill ed circles show the same for apertur es centered on 
CO pe aks. Each point is annotated by its aperture size. A dashed line denotes a best fit to M5 1 by Kennicutt et al. 1 2007). Adapted from Figure 4 of Schruba et al. 
1 2010). Bottom: the simulation equivalent of the top panel in our MC-RTF run at 13.3 Myr into high-resolution evolution. Average locations of local maxima 
found in Figure[8]are displayed. Each point is annotated by the resolution of the performed mock observation. 



size, i.e. higher resolution. This is because at high spatial 
resolution, an aperture at a Ha peak tends not to enclose the 
star-forming regions in cold dense clumps, whereas an aper- 
ture at a H2 peak tends not to include the star-forming regions 
in hot diffuse environment (more to be discussed in £]5.3I >. 
Consequently, Ha emission may not be the best star forma- 
tion indicator at resolutions of < 75 pc. It also implies that 
the spatially-resolved star formation relation in local galaxies 
could potentially be influenced by observation strategies (e.g. 
which peaks to choose, which aperture sizes to use). 

To summarize, we have shown in this section that the 
spatially-resolved Kennicutt-Schmidt relation observed at ~ 
kpc scales breaks down at small scales comparable to the size 
of GMCs, < 75 pc. The correlation varies systematically 
depending on whether apertures are centered on Ha or H2 
peaks. It is because an aperture of a GMC size captures only 
particular stages of GMC evolution. A star formation rela- 
tion of the Kennicutt-Schmidt type is very useful in depict- 
ing the globally-averaged star formation, but may not be as 
effective in describing every single spatially-resolved star for- 
mation events in different clumps and galaxies. A corollary 
to this statement is that one cannot correctly simulate high- 
resolution galaxies if the star formation is simply modeled 
with the Kennicutt-Schmidt relation observed at ~ kpc scales. 

5.3. Evolving Environment Around SFMC Particles 

In order to better understand the various stages and forms 
of star-forming regions, we now turn our attention to the envi- 
ronment around SFMC particles. In Figure[l0] four different 
physical properties in the spheres of 8.4 pc diameter 14 cen- 

14 Or the smallest sphere that encloses the cell a particle belongs to. It 
is almost always a sphere of 8.4 pc except the occasional cases in which 
particles have migrated to a coarsely refined grid. 



tered on all active SFMC particles are drawn as functions of 
particle age, < T < 6 Myr (see Sj2). They include: proton 
number density, gas temperature, proton number density in 
H2, and SFR density from mock Ha emission. The data is 
from a snapshot at 13.3 Myr into the high-resolution evolu- 
tion. For H2 density, the value of which is occasionally neg- 
ligible in some cells, a lower bound of 1 cm~ 3 is imposed for 
a numerical purpose. A red dashed line represents a best fit in 
each log-normal plane. 

In these panels, one may immediately notice the diverse en- 
vironment around SFMC particles, not only throughout their 
lifetime, but even at a fixed epoch of T = Myr. That is, even 
at their birth, the (H2) densities around SFMC particles man- 
ifest a large scatter up to an order of magnitude. It is because 
the stellar radiation from an earlier generation of SFMC par- 
ticles may heat up the star-forming gas clumps, and keep hot 
dense gas > 10 3 cm~ 3 still not turned into stars. The scatter 
in (H2) densities grow larger with particle age, reaching three 
orders of magnitude at T = 6 Myr. One may even observe 
mild evolutions of these quantities as functions of particle age, 
namely: decreasing (H2) density and increasing temperature 
in time. Very young SFMC particles may likely be still en- 
shrouded in cold dense clumps in which they were born; how- 
ever, SFMC par ticles could soon sta rt to drift away from their 
birthplaces (e.g. lOnodera et al1l2010l) . This causes the density 
around a SFMC particle to decrease, and the temperature to 
rise as the particle ages. The trend is fueled by the ionizing 
radiation feedback (0 < T < 6 Myr) and supernova feedback 
(4 < T < 6 Myr) which disperse the surrounding molecular 
gas. Readers should note that this observation is in line with 
the finding in Paper I that the escape fraction from a SFMC 
particle increases on average from 0.27% at its birth to 2.1% 
at the end of a SFMC particle lifetime, 6 Myrs. For these rea- 
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FIG. 10. — Gas density, temperature, H2 density, and SFR density from mock Hce emission in the spheres of 8.4 pc diameter centered on active SFMC particles, 
as functions of particle age. For H2 density a lower bound of 1 cm~ 3 is imposed. A red dashed line represents a best fit in each log-normal plane. 



sons, one can claim that the large scatter and/or the evolution 
of physical properties in the vicinity of SFMC particles are 
responsible for the breakdown of the traditional Kennicutt- 
Schmidt type star formation laws at smaller scales. 

To understand how the star-forming regions are sampled 
with larger apertures (> 300 pc), we now change the reso- 
lution of the previous experiment. In Figure [TT] we plot the 
physical quantities in the spheres of 600 pc diameter centered 
on active SFMC particles as functions of T. Due to the larger 
size of the aperture, averaged properties such as (H2) densi- 
ties look less scattered, and more clustered around the average 
values (red dashed lines). Because several SFMC particles 
tend to congregate in a single star-forming clump that is often 
much smaller than 600 pc, they often exhibit nearly identi- 
cal environmental properties regardless of their ages (i.e. data 
points distributed almost horizontally lined up). This then re- 
sults in little or no evolution of the averaged physical prop- 
erties in each panel (i.e. a flat red dashed lines). Therefore, 
traditional star formation laws still hold in observations made 
with apertures that are large enough to average over various 
evolutionary stages of GMCs. 

5.4. Comparison to the Run Without Stellar Radiation 

So far in Sj5] we have exclusively focused on the MC-RTF 
run which includes both stellar radiation and supernova feed- 
back. Our study is made possible only via the consistent in- 
tegration of high spatial resolution, realistic feedback physics 
and chemistry, and post-production process. To demonstrate 



this, we compare our results in £ )5.2l to that of the MC-TF run 
which includes supernova feedback, but not ionizing radia- 
tion. In Figure [HI we have post-processed the MC-TF run in 
the same manner as the MC-RTF run at the identical timestep 
(lower left panel of Figure[8]). Due to the lack of radiation that 
could have ionized hydrogen in the vicinity of SFMC parti- 
cles, Ha emission and SFR density estimated from it are very 
low in the MC-TF run. Obviously, supernova feedback alone 
fails to achieve the reasonable galactic ISM with a sensible 
density of ionized hydrogen around SFMC particles. One can 
thus conclude that post-processing alone is insufficient to pro- 
duce realistic mock observations like Figure [7] or Kennicutt- 
Schmidt plots. Such maps can be properly generated only 
by self-consistently including radiation calculation in a high- 
resolution galaxy formation simulation. 

6. DISCUSSION 

6.1. Summary and Conclusions 

Using a comprehensive high-resolution simulation of a 
dwarf-sized galaxy including a sophisticated model of stellar 
feedback, we have examined the globally-averaged star for- 
mation, and the spatially-resolved star formation relation on 
a simulated galactic disk. Our goal has been to investigate 
what prompts the spatially-resolved star formation relation 
observed in local galaxies, and why it breaks down at small 
scales. The machinery developed in Paper I is ideal for this 
purpose since it self-consistently computes the interaction of 
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FIG. 1 1 . — Same as Figure [TOl but averaged in spheres of 600 pc diameter. Because most individual star-forming regions are much smaller than the sphere of 
600 pc in diameter, several SFMC particles in the same clump show nearly identical properties of their environment regardless of their ages. 



the ionizing stellar radiation with neighboring gas clouds. Our 
major findings are as follows. 

1. Self-regulated Galactic Star Formation: Our new im- 
plementation of stellar feedback includes ionizing ra- 
diation as well as supernova explosions, and we han- 
dle ionizing radiation by solving the radiative transfer 
equation rather than by a subgrid model ($2}. Joined 
with high numerical resolution of 3.8 pc, the realis- 
tic description of stellar feedback helps to self-regulate 
star formation ( -H.lt . Photoheating by stellar radia- 
tion retains hot dense gas above the density threshold 
for SFMC particle creation, 10 3 cm~ 3 , which otherwise 
would have been unstable against Jeans fragmentation 
and deposited into SFMC particles ( cj4.2| i. 

2. Simulated Observation of Star Formation Relation: Be- 
cause we have self-consistently calculated the location 
of ionized gas, we for the first time are able to make 
spatially-resolved mock observations of galactic star 
formation tracers, such as Ha emission. It is also 
feasible to observe how stellar feedback manifests it- 
self i n the correlation between ionized and molecular 
gas (i )5.11 >. The reported mock observations are possi- 
ble only because the following key ingredients are sys- 
tematically integrated: (a) sufficient resolution both in 
space and time, (b) the realistic stellar feedback model 
including ionizing radiation, combined with primor- 
dial chemistry, and (c) the post-production procedure 



to evaluate Ha emission and H2 fraction on a galactic 
disk. The realistic stellar feedback model established in 
this series of papers provides us with an unprecedented 
insight into the physics of galactic star formation. 

3. Spatially-resolved Star Formation Relation: Applying 
our mock observation techniques to the disk in a galac- 
tic halo of 2.3 x lO n M , we find that the spatially- 
resolved star formation correlation between SFR den- 
sity (estimated from mock Ha emission) and H2 den- 
sity shows large scatter, especially at high resolutions of 
< 75 pc. We also reproduce the phenomenon that cor- 
relation of the Kennicutt-Schmidt type observed at ~ 
kpc scales becomes loose at small scales that are com- 
parable to the size of GMCs, and varies systematically 
depending on whether apertures are centered on Ha or 
H2 peaks. This is because an aperture of GMC size only 
captures particular stages of GMC evolution f <35.2l . 

4. Evolving Environment Around Star-forming Particles: 
By examining the evolving environment around SFMC 
particles, we demonstrate that the large scatter and/or 
the evolution of physical properties in the very vicinity 
of SFMC particles are caused by a combination of stars 
drifting from their birthplaces, and molecular clouds 
being dispersed via ionizing radiation and supernova 
feedback ($53). These factors cause the breakdown of 
the traditional star formation laws at small scales. 
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FIG. 12. — The simulated spatially-resolved Kennicutt-Schmidt plot produced with 75 pc resolution at 13.3 Myr into the high-resolution evolution for the 
MC-TF (left, only with supernova feedback) and MC-RTF runs (right, with both supernova and stellar radiation feedback, same as the lower left panel of Figure 
[8). See the caption of Figure[S]for a detailed description of these plots. Due to the lack of ionizing radiation in the MC-TF run, the SFR density estimated from 
Ha emission is very low in comparison with that of the MC-RTF run. 



6.2. Future Work 

Using the sophisticated feedback scheme developed in this 
series of papers, we plan to work on the following fronts: (a) 
More energetic supernova feedback. In the reported study, we 
adopted a relatively small value for supernova feedback en- 
ergy which is only marginally effective in stopping the unim- 
peded collapse of star-forming gas clumps. While it has been 
a choice to further contrast the effect of ionizing stellar radi- 
ation , we now intend to use more energetic supernova energy 
(e.g. iTasker & Brvanl 120081) . or a more accurate deposition 
of energy using a stellar populatio n synthesis model such as 
Starburst 99 dLeifherer et alJfl999h . From a parameter study 
on different feedback channels, we will investigate the relative 
importance of ionized gas pressure and supernova explosions 
as sources of stellar feedback, (b) Drift of young runaway 
stars. We also aim to include the drift of young runaway 
stars (e.g. ICeverino et alj|20lol) which could potentially en- 
hance the escape of ionizing photons and alter the distribution 
of ionized gas around SFMC particles, (c) Galactic evolu- 
tion. Eventually we plan to perform a galactic simulation with 
sophisticated physics for a longer timescale in order to fully 
comprehend the co-evolution of star-forming gas clumps and 
the entire galaxy. As an intermediate step towards simulating 



a full-blown disk galaxy, we first target a smaller system that 
is comparable to the Small Magellanic Cloud. 
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